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A Fermi accelerator is a billiard with oscillating walls. A leaky accelerator interacts with an 
environment of an ideal gas at equilibrium by exchange of particles through a small hole on its 
boundary. Such interaction may heat the gas: we estimate the net energy flow through the hole 
under the assumption that the particles inside the billiard do not collide with each other and remain 
in the accelerator for sufficiently long time. The heat production is found to depend strongly on the 
type of the Fermi accelerator. An ergodic accelerator, i.e. one which has a single ergodic component, 
produces a weaker energy flow than a multi-component accelerator. Specifically, in the ergodic case 
the energy gain is independent of the hole size, whereas in the multi-component case the energy 
flow may be significantly increased by shrinking the hole size. 


I. INTRODUCTION 


The dynamics of a point particle moving within a 
closed region (billiard) with oscillating walls provide a 
mathematical model for studying the phenomenon of 
Fermi acceleration m\- Such systems typically produce 
an increase in the particle’s kinetic energy and much ef¬ 
fort is devoted to quantify this phenomenon. It has been 
shown that collisions with periodically oscillating walls 
of an ergodic chaotic billiard accelerate the particle so 
that on average its energy grows linearly with the num¬ 
ber of collisions and quadratically as a function of time 
13 [71 Ellin]. In particular, this behaviour is observed in a 
periodically oscillating dispersive billiard Elio] and in a 
stadium with oscillating base [8]. Examples of such bil¬ 
liards are shown in Figs. (a) and (d). It was discovered 
in mm that if the ergodicity of the frozen billiard is vi¬ 
olated, i.e. the shape of the billiard is changed in such a 
way that several ergodic components are created during 
a part of the billiard oscillation cycle, then the average 
energy growth is much faster, typically exponential in 
time. A Bunimovich mushroom deformed so as there ex¬ 
ists particle exchange between its integrable and chaotic 
components corresponds to such multi-component, expo¬ 
nential accelerator [12], see Fig. [^. The mushroom is a 
special case of a large class of billiards with mixed phase 
space where chaotic zones coexist with stability islands; 
the exponential character of acceleration at a periodic 
perturbation of such systems was established in [I31I1]. 
The multi-component accelerators can also be created by 
pseudo-integrability min] and by division of the billiard 
configuration space into disjoint pieces [To], see Fig. [^. 
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One of the primary applications of the Fermi accelera¬ 
tion model is in plasma physics where it is used to study 
the heating of charged particles due to electromagnetic 
waves HITS]. In such systems, the electrons absorb en¬ 
ergy from the wave in the plasma sheath, deposit it in 
the plasma bulk and return back to the sheath. Thus, 
the system is not closed and allows for entry and exit of 
particles. Leaky chaotic systems also emerge in numer¬ 
ous physical situations such as chemical reactions, optical 
microcavities mm and hydrodynamic flows (see recent 
review and references therein m)- 

Stationary leaky billiards have been extensively stud¬ 
ied. It was shown that the escape rate through holes in 
the billiard boundary depends sensitively on the hole po¬ 
sition, size and billiard properties [13 11^231 . In this pa¬ 
per we consider time-dependent leaky billiards and pro¬ 
vide estimates for energy gain for the two above men¬ 
tioned classes of the accelerators: ergodic (like in Figs. 
Ef. c) and multi-component (Figs. [^,d). We consider 
the small hole size limit and observe that TV, the aver¬ 
aged number of collisions a particle spends in the leaky 
accelerator, is inversely proportional to the hole size for 
both cases. Then we demonstrate that for the ergodic 
case the averaged energy gain per particle grows linearly 
with N whereas in the multi-component case the aver¬ 
aged energy gain is much larger and is^^moximated by a 
quadratic polynomial in TV, see Eqs. vs. (14),(16) 

and Fig. [^ 


II. MODEL 

Consider an accelerator which interacts with an ideal 
gas by exchange of particles through a small hole (or a 
few small holes) on its boundary. We assume that the 
gas is at equilibrium, i.e. there is a stationary distribu¬ 
tion of the particles’ speed and the distribution of the 
angles at which the particles move in the gas is uniform. 
We assume that the particles move much faster than the 
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Figure 1: (a) Sin ai accelerator; (b) divided Sinai accelerator; (c) mushroom accelerator; (d) stadium accelerator. See simulation 
section (Sec. |III| ) for geometric specifications and dynamical properties. 





Figure 2: The distribution and the average number of collisions, see Eq. §. (a) dependence on the hole size. Here, the initial 
energy is Ein = 9000 for the Sinai and divided Sinai accelerators, and Ein — 1250 for the mushroom and stadium accelerators, 
(b) no dependence on the initial particles’ energy. Here h = 0.0005 for Sinai and divided Sinai accelerators, and h = 0.00033 
for mushroom and stadium accelerators, (c) the exponential distribution of the exit probabilities for the stadium accelerator. 
Here Ein — 1250 and h — 0.00033. 


billiard boundary. Collisions with the moving billiard 
walls change the kinetic energy of the particles inside the 
billiard, and, on average, this may lead to an outgoing 
energy gain. 
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For a fixed kinetic energy Ein = ^ (we assume the 
particles have a unit mass), the number density of parti¬ 
cles entering the billiard per unit of time is proportional 
to hvin^ where h is the size of the hole (the area of the 
hole for the three-dimensional case). Thus, the incoming 
energy flow at energy Ein is proportional to hvinEin. We 
assume that inside the billiard the particles do not col¬ 
lide/interact with each other, so we can consider each of 
them separately. This gives us the net energy production 
by the accelerator per unit of time: 

Er(^Ein^ — hVin [F/out F/in] 5 (l) 


where Eout is the averaged value of the kinetic energy at 
the moment of exit for a particle that enters the accel¬ 
erator with the energy Ein (we average over all possible 
initial angles and positions in the hole, as well as over the 
phase of the billiard oscillations at the entry moment). 

Let Pj^ denote the probability to exit the accelerator 
after N collisions with the billiards walls and E{N) be 
the corresponding averaged exit energy. Then 

Eout = E Ein)p^ ( 2 ) 

We assume that the hole size h is small enough, so the 
effect of the hole on the statistics of the billiard is negli¬ 
gibly small (as in the case of the stationary Lorentz gas 
[22]b Specifically, we assume that E{N;Ein) can be ap¬ 
proximated by the averaged energy of a particle in the 
closed (i.e. non-leaky) accelerator after N collisions. Ad- 
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Figure 3: Averaged energy gain dependence on the hole 
size. The energy gain increases linearly with Xjh for multi- 
component billiards, see Eq. ( p^ , and is independent of the 
hole size for the ergodic billiards (for sufficiently small holes), 
see Eq. 0. Here, the initial energy is Ein — 9000 for the 
Sinai and divided Sinai accelerators, and Ein — 1250 for the 
mushroom and stadium accelerators. 



Figure 4: Dependence of the average energy gain on the initial 
speed: linear growth with Vin for multi-component billiards, 
Eq.[^ and slow growth with vin for the ergodic case, possibly 
due to order h corrections to Eq.|^ Here h = 0.0005 for Sinai 
and divided Sinai accelerators, and h = 0.00033 for mushroom 
and stadium accelerators. 


ditionally, we assume that and, thus, TV, the averaged 
number of collisions before exit, do not dependent on Ei^ 
nor on the wall velocity u. This is obviously true when 
the billiard walls are stationary, so we extrapolate this 
claim to the case of slowly moving boundaries. We con¬ 
firm this claim numerically for the examples we consider 
here (see Fig. §. In fact, the numerics show that 
can be well approximated by the geometric distribution 
Pn — ~ (1/TV))^“^ (Fig. [^). The average value 

N in this setting is just a geometric characteristic of the 
billiard and the hole. The natural assumption is 
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Figure 5: The speed bias vs. initial speed. Here h = 0.0005 
for Sinai and divided Sinai, and h = 0.00033 for mushroom 
and stadium. 


where h is the size of the hole in the billiard boundary, 
S is the size of the entire billiard boundary, V is the vol¬ 
ume occupied by the billiard, and L is the characteristic 
diameter of the billiard. Relations (§ are confirmed by 
numerical experiments (see Fig. If)- „ 

To find the net energy production Q, it remains to 
estimate the dependence of the averaged energy {E) of a 
particle in the closed accelerator on the number of colli¬ 
sions N. Let us recall how energy is gained in the accel¬ 
erators. The reflection law for a particle hitting a mov¬ 
ing wall is obtained by going to a coordinate frame that 
moves with the same velocity as the wall at the moment 
of collision. In the moving coordinates, we have an elas¬ 
tic reflection law which, after returning to the stationary 
frame, results in the reflection law 

v'^=2u{t,x)-v^, (4) 

where u(t^x) is the normal velocity of the wall at the 
collision point x at the moment t; v and v' are the ve¬ 
locities before and after the collision, and the subscripts 
T and || stand for the components of the velocity which 
are normal and, resp., parallel to the wall. 

If the billiard is chaotic, then the correlations between 
the consecutive angles (j) at which the particle hits the 
wall decay fast. Therefore, the process described by Eq. 

may be approximated by a random walk with reflec¬ 
tions: at each collision the particle velocity undergoes a 
reflection and acquires an increment at a random direc¬ 
tion. 

This random walk proceeds differently for the two main 
classes of accelerators, ergodic and multi-component [7]. 
In the ergodic case the random walk becomes unbiased 
in the large speed limit, which means that the square 
of velocity (i.e. the kinetic energy) grows linearly with 
the number of collisions. Indeed, by taking the square of 
Eq. Q, the energy E^^ = 1'^^ after the N-th collision 
satisfies 

= E,, - 2u{t^,x^)v^ coscp^ +2u^{t^,Xj,). (5) 
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Since \u\ <C the change in the billiard shape and rel¬ 
ative change in the energy are not significant for a large 
number of consecutive collisions, so one can average Eq. 

over the ergodic measure in the (x, 0)-space. The sec¬ 
ond term in Eq. can be much larger than the third 
one, but one can check (see e.g. [7]) that after averaging 
over the ergodic measure and over the period of the bil¬ 
liard oscillation, the second term vanishes (this is a con¬ 
sequence of the existence of the so-called Anosov-Kasuga 
adiabatic invariant in the ergodic case, see miiiiiniiiii- 
EHl). Corrections to the averaging due to a slow change 
in the billiard shape and energy were computed in [2] . It 
follows from [2] that, after the averaging, the 0(?ii;)-term 
in Eq. ^ eff ectively acquires a small factor of order \u\/v 
(see also |lQjb Thus, the effective change of the averaged 
energy per collision is of order (i.e. it is a certain 
portion of the kinetic energy of the wall). 

It follows that in the ergodic case the averaged energy 
of a particle grows as {E{N)) — Ei^ = where u 

is the averaged wall speed (average of |i4|). Thus, we 
conclude (see Eqs. ([^,(|^) that in the small hole limit 

- V 

Eout-Ein = k-N = k — - ( 6 ) 

for some coefficient k that may depend on the billiard 
shape and on the details of the protocol of the billiard 
wall oscillations. By plugging this result into Eq. 
we obtain that the energy gain rate in the ergodic case is 
positive, independent of the hole size h, and is given by 

V 

G{Ein) = kvin--, (7) 


i.e. it is proportional to the kinetic energy of the billiard 
wall, to the volume of the billiard, and inverse propor¬ 
tional to the time Ljvin the gas particle with the speed 
Vin needs to traverse the billiard once. 

Next, we investigate the case of a multi-component ac¬ 
celerator. In this case the ergodicity of the fast motion 
is broken, so the 0{uv) term in Eq. (§ does not average 
out. This means that the random walk in the veloc¬ 
ity space acquires a non-vanishing bias, so the particle 
speed is linear in N and its energy is quadratic in A". A 
more precise description of this process is done based on 
the theory developed in mun]. We note that the time 
between two consecutive collisions tends to zero as the 

L/v^ and it follows 

- Lu^ cos(/)^E;^, i.e. 


particle speed grows, ~ ^ 

from Eq. (5) that ^ 

|_J ^iV + l 

in the non-ergodic case the energy changes exponentially 
with time, with a certain random rate. On a longer time 
scale this process can be modeled by a multiplicative ran¬ 
dom walk (see [71 [TOl[T4] b 
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where E^ = '^ is the kinetic energy after n periods of 
the billiard oscillations, and is the sequence of inde¬ 
pendent, identically distributed random variables, inde¬ 
pendent of the initial energy. Importantly, it is shown 


in mill] that this random walk cannot be decelerating 
and, typically, 

Elne„>0, Ee„>l, E^2^1. (9) 

If we ignore details of particle behavior on the time 
scales below the period T of billiard oscillations, we can 
infer from Eq. ^ the following description for the be¬ 
havior of the averaged energy and speed gain at time t: 

{E{t)) = EinC^^, {v{t)) = (10) 


(^^(S1)^^(S2)) = {v\s,)) ■ {v{S 2 )/v{s,)) = 

= (s2 > si), 


where /i = ^ InE^^ > 2A = InE^ > 0. 

Note that the number of collisions up to time t can be 
related to the particle speed via LN(t) ^ v(s)ds. So, 
the averaged number of collisions up to time t is given by 


(N(t)} = kiVin 


(e^^ - 1) 

XL 


(see Eq. diol )) and, by Eq. ( |IT| ), 

{N^{t)) = f f {v{si)v{s2))dsids2 
Jo Jo 


h 

£2 

£2 


n {v{si)v{s2))ds2dsi 
1 


= 4fc2 


Ei, 


(m - A)L2 


e#** -1 - 1 


where ki ^2 are some coefficients of order 1. 

Rearranging the expressions for {N) and {N'^), we get 


{E{t))-Ein = 2kiixVinL{N {t))+Ak2lJ? 



{Nitf), 


which, after averaging over the time t the particle resides 
in the billiard, gives 


_ L2 

Eout Ein ^ 


2ki^^N+Ak2{l- ^ 


( 12 ) 


The rate fi in this formula is a well defined quantity de¬ 
termined by a one period run of the accelerator, and an 
analytic expression for /i is also available in many cases 
[aiMi]. However, to compare the energy gain with that 
given by Eq. in the ergodic case, we need to relate 
the rate /i with the kinetic energy of the wall. 

ByEq. (§ the exponential growth rate /i is always 
non-negative. The minimal value /i = 0 is achieved when 
the distribution of particles in the billiard remains uni¬ 
form during the period of billiard oscillations, like in the 
case of ergodic billiard. Hence, when the deviation from 
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the ergodic behavior is small, the rate /i is small and be¬ 
haves like the square of a certain quantitative measure 
of this deviation. The violation of ergodicity in the ex¬ 
ponential accelerator is caused by changes of the phase 
space structure of the frozen billiard as its shape changes 
with time. Therefore, we relate the deviation from er¬ 
godicity to the magnitude of the shape change over the 
period. As the billiard size is changed with the mean 
speed tZ, the dimensionless parameter estimating the er¬ 
godicity violation is uTjL. Thus, the rate of the energy 
increase over the period T is, at small tZ, given by 

(13) 


One can also extract this relation from the formulas for 
/i for various cases of multi-component accelerators (e.g. 
from m H uni na US]). Notice that even when u is not 
small, the ratio uTjL remains bounded and so does {iT 
(we assume everywhere that the typical length scale of 
the billiard does not change hugely along the cycle; oth¬ 
erwise pathological behaviors may arise). So, Eq. (13) 
can be used in this case as well; it would then simply 
mean that the quantities on both side of the relation are 
of order 1. 


Plugging Eq. (13) into Eq. (12), we find 


■72 r 


A'out ^in — 


ki 


+ k2iiTN‘^ 


(14) 


or, in the case uTjL 1, 

VinT 


V? 

^out E/in ~ 


+ iV2 


, (15) 


where the new coefficients ki ^2 depend on the shape of 
the billiard and the protocol of the wall oscillation. 

The term VinTjL is proportional to the number of col¬ 
lisions per period. In our setting this number is assumed 
to be large. It is also important for the validity of the ex¬ 
ponential growth model that the particle is initially fast 
(i.e. Vin ^ u) yet it remains in the billiard at least for 
one period. Therefore, Eqs. (12), (14), and (15) are valid 
under the assumption 


l«%^< AT; 


in particular < ^. When this condition is violated, 
our theory is not applicable. 

Comparing Eqs. ( p^ , (15) with Eq. ^ we see that the 
energy gain in the multi-component accelerator is much 
larger than in the ergodic case. Even if the exponential 
growth rate /i is very small, the coefficient of u‘^N in 
Eq. (is large whereas the corresponding coefficient in 
Eq. (^ is simply a constant. With the increase of the 
quadratic in N term becomes dominant in Eqs. (14), (15) 


and provides the main contribution to the energy gain. 


Using Eq. (§ for N and we finally find the energy 
production rate of Eq. 0 for the multi-component case: 


g(e;,o = 


0-72 


'^in i- 


Vu 


L 2 



+ k2iiT 


V_ 

Lh 


(16) 


Clearly, the gain rate G can be made much larger than 
the gain in the ergodic case by diminishing the hole size 
or by increasing the incoming velocity. 


III. SIMULATIONS 

There are two distinct predicted dependencies of the 
energy gain on the hole size and on Vin for the two types 
of leaky accelerators. To examine these predictions we 
consider two classes of leaky accelerators where each of 
these is considered with two sets of parameters - one cor¬ 
responding to an ergodic case and the other to a multi- 
component case. 

Dispersing accelerators (Eigs. 0^ ,b): at t = 0 a 
vertical bar is inserted at a position to the double 
Sinai billiard (a rectangle with two discs). Then the bar 
moves to the right with a constant velocity u till time 
r and then the bar is removed. The cycle restarts at 
time T. We consider two cases: la) ergodic case where 
the bar only partially blocks the rectangle, covering 90% 
of the rectangle length (this case is called “Sinai” in all 
Eigures); lb) multi-component case where the bar com¬ 
pletely divides the rectangle into two parts (hereafter “di¬ 
vided Sinai”). Notice that each component of the frozen 
billiards is ergodic and mixing m- These accelerators 
exhibit exponential-in-time energy growth in the multi- 
component case and quadratic-in-time energy growth in 
the ergodic case m- To examine the leaky behavior, 
two holes of length h are placed on the upper rectangle 
boundary; the holes are shifted from the disks centers 
to avoid fast escaping orbits, and two holes are intro¬ 
duced to avoid strong dependence on the billiard oscil¬ 
lation phase. In all simulations we use the following pa¬ 
rameters: the rectangle width is a = 4, its height is 6 = 2, 
the disks radii are 1/2, the bar velocity isu = 0.1, the bar 
is introduced at the position = 0.0915 and removed 
at the time moment r = 1.83, the period is T = 5.49. 
The initial energy and the hole size are as indicated in 
the figures. 

Focusing accelerators (Eigs. [^,d): The mushroom 
is a multi-component system having an integrable com¬ 
ponent and a chaotic component [29] whereas the slanted 
stadium is ergodic and mixing m- The oscillating 
mushroom accelerator exhibits exponential-in-time en¬ 
ergy growth [T2| whereas the oscillating stadium exhibits 
quadratic-in-time energy growth laiHi. 

The shape of the mushroom is determined by the fol¬ 
lowing four parameters: r is the radius of the cap; w is 
the half-width of the hole at the bottom of the cap, it 
coincides with the half-width of the stem at its highest 
point and w < r; £ is the length of the stem; and the 
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angle 0 describes the inclination of the stem sides. When 
w = r the mushroom becomes a slanted stadium. 

For the purpose of numerical experiments we used the 
following protocols: r{t) = 1, w{t) = bo — bi{l — cos(t)), 
and i{t) = ao — ai sin(t), uq = 1, 6o = 1. For the mush¬ 
room, we set ai = 0.5 and bi =0.4; and for the stadium 
ai = 0.5, bi = 0. The hole is located at the bottom of the 
stem with the center displaced by 0.01 from the center of 
the stem. In all experiments 0 = 0.1111. 

In each numerical experiment 2000 particles are in¬ 
jected at randomly chosen times during the period [0, T] 
into the billiard through the holes, with random posi¬ 
tion in the hole and entering angle. Each particle moves 
inside the billiard undergoing elastic collisions with the 
boundary till it exits by colliding with the hole. The exit 
time, the number of collisions till exit, and the exit speed 
are recorded. 

Figure demonstrates that the average number of 
collisions at exit, scales linearly with l/h for both the 
multi-component and the ergodic cases. Figureshows 
that N does not depend on the initial energy, again for 
both cases. Figure shows the geometric distribution 
of the exit collisions for the mushroom. These results 
support the assumptions made in Eq. <§• 

Eigure shows the dependence of Gjvin = (Eout — 
Ein)h onl/h for the four billiard types. Eigure^ shows 
that the net energy flow increases linearly with l/h for 
billiards with exponential acceleration (non-ergodic case) 
as predicted in Eq. (16). The inset shows that for suf¬ 
ficiently small holes the flow is essentially independent 
of the hole size for the ergodic billiards as predicted by 
Eq. 0. 

Figure 1^ shows that the average energy gain, Eout — 
Ein^ grows linearly with Vin for the multi-component 
cases as predicted by Eq. (14). The inset shows that 
this difference also grows with Vin for the ergodic case 
(at a slow pace). To the leading order in h one should 
not see such growth according to Eq. §• We explain 
the effect by the order h difference between the statistics 
of the closed and leaky accelerators. Eor example, aver¬ 
aging the second term of Eq. (§ over a boundary with 
the order h hole produces corrections of order huv^ which 
cause an order h bias in the random walk of the velocity. 
Eigure 1 ^ shows the mean bias {vout — in) / {N)^ which is 
indeed present and does not vanish in the limit of large 
initial speed, both in the ergodic and multi-component 


case. However, the bias is much smaller in the ergodic 
case, confirming our conjecture that it is of order h. The 
dependence of the energy gain on u appears to be a more 
delicate numerical issue and will be further investigated 
elsewhere. 


IV. DISCUSSION AND CONCLUSION 


A priory, a faster accelerator does not automatically 
imply higher energy gain as faster particles escape earlier 
from the accelerator and the net energy gain reflects the 
balance between the averaged escape time and the energy 
gained by then. In this paper we find that in a billiard 
the escape is mainly determined by the average number 
of collisions, and from this observation we are able to find 
the balance between these two factors for two different 
types of accelerators. Our results conclusively show that 
multi-component leaky accelerators produce higher net 
energy flow than their ergodic counterparts. Moreover, 
we show that the difference increases significantly when 
the hole size decreases. 

One may worry that decreasing the hole size will lead 
to very long residence times. However, for the multi- 
component case the residence time increases only loga¬ 
rithmically with the hole size. Indeed, the number of 
collisions is related to the residence time via LN(t) ^ 
Jq v{s)ds^ hence, in the multi-component case the av¬ 
eraged residence time is proportional to log TV logh, 
while in the ergodic case it is proportional to Vn ^ 
^ Hence, we expect that taking the small hole limit 
for gaining energy may become practical in the case of 
multi-component leaky accelerators: they produce high 
gain at a fast pace. 
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